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Abstract 

We propose a model for neutrino mass generation in wich no physics beyond a TeV 
is required. We extend the standard model by adding two charged singlet fields with 
lepton number two. Dirac neutrino masses m UD < MeV are generated at the one loop 
level. Small left handed majorana neutrino masses can be generated via the seesaw 
mechanism with right handed neutrino masses Mr are of order TeV scale. 
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Solar and atmospheric neutrino oscillation experiments 0, give a strong evidence for 
the neutrino masses and mixings. In the minimal standard model, neutrinos are massless, 
due to the abscence of the right handed neutrino and the conservation of lepton number. 
Non zero Dirac neutrino masses can be obtained by adding, singlet fermions (right handed 
neutrinos) to the SM. However, in this case one needs to fine tune the Yukawa couplings to 
be ~ 10~ 12 , which is quite unnatural. 

An elegant and natural explanation of the smallness of the neutrino masses compared 
to the quarks and charged lepton masses is provided by the seesaw mechanism ||. In this 
case the Majorana neutrino mass matrix is given by: 

m u = -m^M^ l m D , (1) 

where is the Dirac neutrino mass matrix and Mr denotes the right handed Majorana 
neutrino mass matrix. This mechanism can be implemented naturally in many extensions 
of the standard model gauge group (e.g, in SO (10) the right handed neutrino masses are 
associated with the (B — L) symmetry breaking scale) JJ. The smallness of the neutrino 
mass is due to the suppression by the scale of the heavy fields (right handed neutrinos ) 
after integrating them out. An equally attractive mechanism has been proposed recently 
|H where the lepton number is broken expilicitly by the coupling of a heavy higgs triplet to 
the standard model doublet. 

Another interesting way to generate small Majorana neutrino masses is given by the 
Zee model||, in which the masses are generated at the one loop level. In addition to the 
standard model higgs doublet §sm the model contains a charged scalar field h!- + \ that is a 
singlet of SU(2)l, and another doublet scalar field $ which has the same quantum numbers 
as the standard model higgs. The singlet carries lepton number —2, so that the total 
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lepton number L is conserved in the Yukawa sector. The lagrangian of the Zee model reads: 

Czee = UL T a ClT 2 Lph^ + ^ T SM lT 2 §h^ + h.C, (2) 

where L a is the lepton doublet, C denotes the charge conjugation, and ir 2 is the SU(2) 
Levi-Civita symbol. Lepton number is broken explicitly by two unites throught the cubic 
coupling fx^ M ir 2 ^h^\ In addition, a discrete symmetry is needed such that the higgs 
doublet $ does not couple to leptons. A calculable neutrino Majorana masses are generated 
at one loop through the exchange of the physical higgs and charged lepton fields. Fermi 
statistics imply that the coupling matrix f a p is antisymmetric, which leads to a neutrino 
mass matrix with vanishing diagonal elements. 

The possibility of applying the radiative mechanism to generate Dirac fermion masses, 
including neutrinos, have been studied previously in the context of left right symmetric 
and superstring inspired models ||. The authors of Ref. J7| have presented a model with 
large magnetic moment of the electron neutrino, by adding a singly charged higgs field , in 
addition to the right handed neutrinos. However, in this model the one loop neutrino Dirac 
mass is logarithmically divergent and must be cancelled by adding a counterterm, making 
the neutrino masses arbitrary. 

In this letter we propose a model in which lepton number is broken by the right handed 
Majorana mass Mr ~ TeV, and where calculable Dirac neutrino masses are induced at one 
loop via the exchange of singlet charged scalar fields. 

The model is the standard model extended by introducing two singlet charged scalar 
fields, Si and S 2 . No extra doublet is needed here. We also introduce three families of 
right-handed neutrinos z/r x , z/r 2 , z/r 3 . Let us assume for the moment that the total lepton 
number is conserved. The most general relevant part of the lagrangian consistent with this 
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symmetry is given by: 

C ext = f a pLlCiT 2 LpS{ +) + g a pl T R CvRA +) + h - c ( 3 ) 

where, a, f3 denote generation indices. The couplings f a p are antisymmetric in a and (3, 
while g a p are arbitrary. Conservation of lepton number requires that s[ + ^ and S 2 + ^ carry 
lepton number L = — 2. However at this stage a coupling between the standard model higgs 
and the right handed neutrino is allowed. The presence of such term in the lagrangian leads 
to a Dirac mass for the neutrino, and one needs to fine tune the Yukawa coupling to produce 
m v <C m e . To forbid these couplings we impose a Z 2 symmetry acting as follows: 

(L a , $ S m, S{ +) ) — (L a , $ 5M , s[ +) ) (4) 

(N Ra ,S { 2 +) )^-(N Ra ,si +) ) (5) 

In the limit where Z 2 is unbroken the Dirac mass for the neutrino vanishes to any order of 
perturbation theory. We will therfore assume that Z 2 is broken softly in the higgs sector. 

The scalar potential is defined by: 

V = V^sm&sm, SiS\, S 2 Sl) + kS^I (6) 

where Fi($sm^smi SiSi) contains the usual terms with positive mass square and quadratic 
terms for Si and S 2 such that < Si >=< S 2 >= (since they carry electric charge). The 
term kSiS 2 is the Z 2 soft breaking term which can be of order of the charged singlet higgs 
masses m s . The physical charged higgs fields are given in the basis S[ + ^ and S 2 + ^ by: 

hi = s[ +) cos6 + S { 2 +) sine 
h 2 = -S ( i +) sin6 + Si +) cos6 (7) 

where 
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Dirac masses of the neutrino are generated from the one loop diagrams shown in Fig.l and 
in the basis in which the charged lepton mass matrix is diagonal are given by 

m v a p = ga-ym-ffjp sin 6 cos 9I(M 1 , M 2 , mp) (9) 

where mp are the masses of the charged leptons, Mi, M 2 are the masses of the physical 
charged higgs hi and /i 2 and /(Mi, M2, mp) is the one loop integral given by: 

in evaluating this integral the charged lepton masses can be safely neglected so that the one 
loop induced Dirac neutrino masses are given by 



1 M 2 2 

m v af3 = -^9*i m ~tfiP sin ( 2 0) ln (^) ( n ) 



In this model neutrino masses are proportional to the charged lepton masses, while in the Zee 
model the Majorana neutrinos are quadratic in the charged lepton masses. Note that one of 
eigenvalues of the neutrino mass matrix in Eq. (11) vanishes, since the 3x3 antisymmetric 
matrix f a p has zero determinant. The same feature has been noticed in another model || 
with a doubly charged scalar. The higher loop correction can always be written in the form 
g T Yf, where Y depends on the parameters of the model. Thus at any higher loop level one 
of the neutrinos remains massless. Note that although the charged lepton masses receive 
corrections at one loop by replacing the charged lepton fields in Fig.l by the Dirac neutrino 
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fields, this correction is 5rrii ~ which can safely be neglected. 

Now let us study the phenomenological constraints on the parameters of the model. The 
physical charged scalar fields hi and hi can mediate the decays fi — > ev e v^ r — > ev e is T , 
and r — > /-tz^iv The muon decay (see Fig. 2) can be described by the effective lagrangian: 

£eff = yl"^ 1 + S L^)e L J P ^e L ^^Y^L + $ R^rY V eR V m Y fl R + 5 LR ^ C R V eL V m \Lj\ , (12) 



where 



f 2 
All 

V2G F M Z 



<w = J:» 2 , (is) 



,2 



= (14) 

^ = M?< (15) 



with: 

f2 _ M 2 M 2 

Ml cos 2 + Mf sin 2 



For the r — > ez/ e z/ T and r — > [iv^v^. decays one has to substitute 5 efl by 5 er and 5^ T 
respectively The more stringent constraint on f efl is derived from universality of the \i 
decay which gives @ 

4r < lO^GeV' 1 . (17) 

M ~ K J 

A useful constraint on g efl can be obtained from e — [i universality which gives the bound 

9 -^k E < 3 x 10- 6 GeV- 2 . (18) 
M 

The charged scalars hi and hi contribute to the /i — ► e^y . A straightforward calculation 
gives: 

1 vrc* 

rGu — ery) ~ | (//+) 12 + (^ + )i2 | 2 (19) 

3847T 4 

and by using the experimental upper limit B(fi — > ey) < 1.2 x 1CT 11 we obtain 

I (// + )i2 + (^ + )i2 |< 2.8 x 10- d GeV- 2 . (20) 

M 

Another bound on the (V + A) couplings g a p can be drived from the constraint on the 



Michel parameter £. Using the experimental bounds from Ref. |T(| , we obtain 

a 2 

^<3x l(T 4 GeV- 2 (21) 
M 
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Note that the decay fi — > V e v^ R eR will be forbidden kinimatically in the seesaw case, since 
the mass of the right handed neutrino is much bigger than the mass of the muon. Note that 
for Mi 2 — lOTeV and / < 1, g < 1 the above constraints are satisfied. 

Now let us turn to the mixings and the masses relevant for neutrino oscillations, we 
assume for simplicity g a p ~ f a p, such that the model has four remaining parameters. In 
this case the neutrino mass matrix is symmetric, and has the following form: 



m,. 



m 



f 2 + f 2 

J er 1 Jefj,\m T 
fer f/ir 



f erf /it fefif) 
far feufe 



lefiJeT\ mT J 



(22) 



_f f ( Trigs f r (m^\ f2 ( rn R \ 

where uiq corresponds to the scale of neutrino mass in this model, and it is given by: 



An 2 



rriQ 



-m T - 



(23) 



64tt 2 t Mf + 4k 2 ' 

where in deriving the above relation we have assumed that the charged higgs h% is heavier 
than h\. For k ~ M 2 , m ~ 2Mey, which implies that the parameters f a p < 10 -3 . It is well 



known that a stable MeV neutrino mass can be problematic for cosmology [11 1. Morever 



primordial nucleosynthesis excludes a Dirac neutrino mass from .3 to 25Me^ |fL2|| , which 
implies that f a p should be less than one. 

If there is no strong hierarchy between the Yukawa couplings f a p, then the matrix 
element be safely neglected, and the eigenvalues of m v are 



m 



m U2 ~ 



m U3 ~ 







1 ,m,. 



(24) 



This means that m V3 is the scale relevant for atmospheric neutrino oscillations, which implies 
that ~ f eT ~ 10~ 4 . However the explanation of the solar neutrino anomaly requires a fine 
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tuning between f eT and f^ T such that m U2 << Am 2 tm . The matrix m u can be diagonalized 
(neglecting CP violation in the leptonic sector) by the following orthogonal transformation: 



V 



2ad 
1 



a+b+Va 2 +b 2 
a— 6— \/d z +b' 2 

1 

C 2 



a+b-Va 2 +b 2 
a-b+Va 2 +b 2 

1 

Cz 



(25) 



where 



a = 

b = 

Cl = 

c 3 ~ 



JilT 



e fe/i 

VI + 4a 2 



2(a 2 + b 2 + aV^Ttf + W 2 

4 

2(a 2 + b 2 + bVti r Tb 2 ~+-)^ 

4 



(26) 
(27) 

(28) 
(29) 
(30) 



From the above mixing matrix one can see that the large mixing angle required for the atmo- 
spheric neutrino oscillations can not be satisfied simultanously with the required smallness 
of the U e 3 element. One could ask if the model can accomodate both the atmospheric and 
solar neutrino anomalies for g a p ^ f a p. However we have performed an analysis of all pos- 
sible parameters f a p, g a/ 3, and we found no solution that can fit both atmospheric and solar 
neutrino puzzles with [7 e3 < 0.2. 

Now let us assume that lepton number is broken by the right handed neutrino mass 
term: 

(31) 



£al=2 = ^M Rap N T R ^CN R ^ + h.c 



The masses of the light neutrinos are obtained by diagonalizing the following mass matrix 
in the basis vl, N r : 

rn M°A < 32 > 
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where m UD is the Dirac neutrino mass matrix induced at one loop as in Eq.(ll). The seesaw 
approximation (Mr » m UD ) implies 

M (n 9 ht) „ _ m T D M- l m UD (33) 

We will assume that the Yukawa coupling constants f a /3, g a /3 are of order unity, and Mr, Mi 
and M 2 in the TeV range. In this case the light Majorana neutrino masses are smaller than 
eV. Of course, one of the eigenvalues will vanishe, and this requires the other masses to be 
in the atmospheric and solar neutrino ranges. This can be easily accomodated since the 
seesaw fomula is very sensitive to any deviation of the coupling constants from unity, since 
it varies as the fourth power of the coupling constants. 

The model can not predict the mixing angles since the parameters f a p and g a p are 
arbitrary, in contrast to the Zee model, in which the neutrino mass matrix has vanishing 
diagonal elements without any fine tuning or imposing extra symmetries. However, unlike 
the Dirac case, the seesaw neutrino mass matrix can accomodate the atmospheric and the 
solar anomalies. Choosing the Yukawa parameters j a p^g a p to take the values: 

9n = #22 = g-33 = 92i = gi3 = 531 = 1 

#23 = #32 = 

g 12 ~ -.23 

/23 = 1 (34) 
f 12 = -f 13 ~ .5 (35) 

with the right handed neutrino masses (Mr 1} Mr 2 , Mr 3 ) ~ (15.2, 3.3, 10)TeV, leads to a 
bimaximal mixing with large angle MSW [|HJ effect as the solution to the solar neutrino 
puzzle. The numbers above suggest that one can easily choose a texture for the matrix g a p, 
that can fit both the atmospheric and solar neutrino data. 
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In summary we have proposed a model for generating small majorana neutrino masses 
at the TeV scale using the seesaw mechanism. The smallness of the neutrino masses in this 
model is not due to the largness of the right handed neutrino mass, but to the smallness of 



the Dirac masses. The model may also have an intersting implication for cosmology |L4 
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Figure 1: One loop diagram giving rise to Dirac neutrino masses 



11 




12 



